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I.  FORWARD 


The  Workshop  on  Constitutive  Modeling  of  Air  and  Water  Saturated  Sand  for  Shock  Prop¬ 
agation  Modeling  was  held  on  September  20-  21,  1994,  at  the  Naval  Research  Laboratory, 
Washington,  DC.  The  workshop  was  sponsored  by  the  Naval  Research  Laboratory  and  was 
chaired  by  Mark  H.  Emery,  Code  6440,  and  Philip  J.  Valent,  Code  7401.  The  objectives  of 
this  workshop  were  to  appraise  the  present  state  of  knowledge  with  respect  to  the  character¬ 
ization  of  the  air  and  water  saturated  sandy  medium,  the  experimental  evidence  for  shock 
mitigation  and  propagation  in  this  medium,  the  present  constitutive  modehng  capabili¬ 
ties,  and  to  develop  a  focused  research  effort  to  address  the  complex  physics  issues  related 
to  this  task.  Twenty-five  scientists/engineers  attended  the  workshop  with  participation 
from  NRL,  NSWC,  ONR,  DNA,  the  university  community,  and  private  industry.  Partici¬ 
pants  discussed  issues  associated  with  mine  cormtermeasme  systems  including  shock  wave 
propagation  in  sandy  sediments  with  water  and  free-air  pore  fluid,  measurement  of  the 
mechanical  properties  of  both  dry  and  partially  satmated  sand,  hydrodynamic  modeling 
of  this  environment,  and  continuum  and  micromechanical  aspects  of  consititutive  models. 
Subsequent  sections  provide  the  workshop  announcement,  planned  agenda,  summaries  of 
the  participants’  presentations,  and  findings  and  recommendations  (incorporating  informa¬ 
tion  shared  in  open  discussions  with  the  participants)  which  includes  the  opinions  of  the 
workshop  chairs  (Emery  and  Valent)  which  may  differ  from  the  opinions  of  the  workshop 
participants. 

This  workshop  serves  only  to  initiate  the  scientific  and  technological  (S&T)  dialogue  ad¬ 
dressing  the  research  and  development  requirements  of  the  Navy  and  Marine  Corp.  Prom 
this  beginning  it  is  anticipated  that  some  focussed  set  of  S&T  programs  will  be  formu¬ 
lated.  These  will  involve  a  broader  spectrmn  of  scientists  and  engineers  whose  expertise  is 
required  to  address  the  relevant  problems.  The  need  for  research  in  this  area  is  compelling 
and  we  seek  to  lead  that  research  in  directions  offering  the  most  promise  for  the  DoD  while 
advancing  the  state-of-the-art  benefiting  the  nation  as  a  whole. 
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II.  WORKSHOP  ANNOUNCEMENT 


WORKSHOP  ON  CONSTITUTIVE  MODELING 
OF  AIR  AND  WATER  SATURATED  SAND 
FOR  SHOCK  WAVE  PROPAGATION  PREDICTION 

September  20  -  22,  1994 

The  US  Naval  Research  Laboratory 
Washington,  DC  20375-5000 

Workshop  Sponsored 
by 

Naval  Research  Laboratory  and  Office  of  Naval  Research 

Background  and  Aim: 

The  US  Navy  is  in  urgent  need  of  a  self-consistent  numerical  capability  for  modeling  explosion 
and  shock  effects  in  surfzone  seafloor  sediments,  primarily  sands,  with  seawater  and  free-air  pore 
fluid.  One  of  the  most  complex  tasks  in  this  effort  is  the  development  of  accurate,  efficient,  and 
physics  based  constitutive  models  for  gassy  sediments  in  and  near  surfzone  environments.  The 
gassy  sediment  environment  is  a  complex  microscale  medium  characterised  by  sediment  structure 
parameters  (porosity;  permeability;  tortuosity;  grain  size,  shape,  mineralogy,  and  packing;  shear 
strength;  anisotropy),  gas  parameters  (bubble  size,  shape,  orientation,  and  distribution;  gas  vol¬ 
ume),  and  fluid  parameters  (distribution;  water  volume;  compressibility).  Improved  constitutive 
models  are  needed  for  progress  in  performance  prediction  of  mine  countermeasure  efforts. 

The  aims  of  this  workshop  are  to  appraise  the  present  state  of  knowledge  with  respect  to  the 
characterization  of  this  medium,  the  experimental  evidence  for  shock  mitigation  and  propagation 
in  this  medium,  the  present  constitutive  modeling  capabilities,  and  to  develop  a  focused  research 
effort  to  address  the  complex  physics  issues  related  to  this  task. 

Proposed  Topics: 

Characterization  of  the  medium  (in  support  of  shock  propagation  modeling) 

Presentations  will  focus  on  detailed  descriptions  of  the  gassy  sediment  environment,  emphasizing 
sand  sediments,  the  degree  of  variability  of  the  sediment,  gas,  and  fluid  parameters,  the  relative 
physical  importance  of  these  variables,  and  the  resolution  of  the  techniques  for  measurement  of  the 
parameters  for  input  for  the  constitutive  models. 

Acoustic  and  shock  wave  propagation 

This  topic  will  focus  on  the  availability,  resolution,  and  accuracy  of  experimental  results  (both 
in  situ  and  laboratory)  pertaining  to  acoustic  and  shock  wave  propagation  in  this  medium.  The 
relative  roles  of  the  physical  variables  (and  their  sensitivity)  that  influence  the  propagation  char¬ 
acteristics  will  be  emphasized. 

Constitutive  modeling 

Emphasis  will  be  placed  on  clarifying  the  present  modeling  capabilities,  the  physics  that  is  presently 
being  modeled,  deficiencies  in  the  present  modeling  capabilities,  and  potential  solutions  to  these 
problems.  Topics  of  interest  include  microscale  modeling,  macroscale  modeling,  and  the  synthesis 
of  macroscale  models  from  microscale  models. 
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III.  AGENDA 


Workshop  on  Constitutive  Modeling 
of  Air  and  Water  Saturated  Sand 
for  Shock  Propagation  Prediction 

September  20-21,  1994 

The  US  Naval  Research  Laboratory 
4555  Overlook  Dr  SW 
Washington,  DC  20375 


Tuesday,  September  20,  1994 


8:30  AM 

9:30  AM 
9:45  AM 
10:15  AM 
10:45  AM 
11:15  AM 
12:00  PM 
1:30  PM 

2:00  PM 

2:30  PM 
3:00  PM 
3:30  PM 

4:00  PM 


Registration  and  Coffee 

M.  Emery  (NRL),  Welcome  and  Opening  Remarks 
H.  Chen  (NSWC/WO),  Overview  of  the  Mine  Counter  Measure  Problem 
R.  Barsoum  (ONR),  Constitutive  Modeling  of  Sands/Soils 
P.  Valent  (NRL/SCC),  Characterization  of  the  Medium 
Discussion  Session 
Lunch 

R.  Dick  (Univ  MD),  Hopkinson  Bar  and  Explosive  Studies  of  Wet  and 
Partially  Saturated  Sands 

D.  Chitty  (ARA,  Inc.),  Field  and  Laboratory  Investigations  of  Wave  Propagation 

in  Saturated  Sand 

R.  Walker  (ARA,  Inc.),  Ground  Shock  Attenuation  in  Three  Phase  Earth  Materials 
Coffee  Break 

R.  Pandey  (Univ.  S.  Miss.),  A  Microscopic  Simulation  Approach  to  Fluid 
Flow  Through  Porous  Media 

Smmnary  and  Discussion 


Wednesday,  September  21,  1994 


8:30  AM 
9:30  AM 

10:00  AM 

10:30  AM 

11:00  AM 

11:30  AM 
12:00  PM 
1:30  PM  - 


Registration  and  Coffee 

A.  Shukla  (Univ  RI),  Experimental  and  Numerical  Studies  of  Wave  Propagation 
in  Granular  Materials  Subjected  to  Explosive  Loadings 
D.  Flippen  (NRL),  Methods  of  Synthesizing  Macroscale  Models  From  Microscale 
Models  for  Composites  with  Planned  Adaptions  to  Sand 
D.  Zhang  (Johns  Hopkins  Univ),  Ensemble  Phase  Averaged  Equations  and 
Constitutive  Relations  for  Dispersive  Flows 
D.  Peric’  (Univ.  Co-Denver),  Discontinuous  Bifurcations  in  Partially  Saturated  Soils 
Under  Undrained  Conditions 

R.  England  (TRT,  Inc),  Review  of  Previous  DNA  Water  Shock  Program 
Limch 

Summary  and  Discussion  (Open  Workshop) 
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In  all  three  topics,  the  focus  of  the  workshop  will  be  on  the  physical  properties  of  this  complex 
medium,  assessment  of  the  present  state  of  the  art  in  modeling  those  properties,  and  potential 
improvements  to  existing  models  or  development  of  new  constitutive  models. 

Format: 

The  format  of  the  workshop  will  consist  of  two  full  days  of  thirty-minute  presentations  (plus 
discussion  time)  by  workshop  participants.  The  agenda  will  be  separated  into  the  three  general 
topics  discussed  above.  There  will  be  a  thirty-minute  discussion  period  at  the  end  of  each  general 
session.  In  addition,  one-half  of  the  third  day  will  be  a  roundtable,  brainstorming  session  for  all 
participants  which  will  review  the  present  state  of  the  art,  assess  the  physics  issues,  their  relevance, 
the  means  to  measure  them  and  model  them,  analyze  the  present  constitutive  modeling  capabilities, 
and  develop  a  research  agenda  (both  experimental  and  numerical)  to  address  the  physics  issues  in 
order  to  accelerate  the  progress  of  performance  prediction  of  mine  countermeasure  efforts. 
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IV.  SUMMARY  OF  PRESENTATIONS 


Chen  and  Barsoum  provided  an  overview  of  the  mine  countermeasure  problem  from  the  6.2 
and  6.1  perspective  respectively.  The  concept  entails  the  projection  of  an  explosive  net  onto 
the  surf  zone  and  subsequent  detonation  of  the  net  to  clear  a  safe  path  for  landing  craft. 
This  procedmre  is  simple  in  concept  but  complex  in  execution.  To  predict  performance  of 
such  a  system  requires  the  ability  to  predict  shock  propagation  through  wet  sands  that 
are  partially  filled  with  gas.  The  performance  of  this  system  is  being  estimated  by  field 
tests  in  test  pond  sites  and  by  computational  simulations  as  there  can  be  no  tests  of 
the  explosive  array  on  real  beaches  based  on  cmrent  environmental  policies.  In  addition, 
the  distributed  explosives  require  high  energy  and  low  charge  diameter.  Chen  presented 
results  from  a  two  dimensional,  -  based  CTH  simulation  of  a  charge  array  on  a  sand 
surface  under  water.  He  also  discussed,  conceptually,  a  smooth  particle  hydro  (SPH) 
model  as  a  micromechanical  approach  to  modehng  the  sand.  Barsoum  discussed  several 
6.1  endeavors  -  split-Hopkinson  bar  experiments,  micromechanically  based  constitutive 
models,  and  double  asymptotic  approximations  for  wet  sand/mine  interaction.  Valent 
discussed  several  aspects  of  the  characterization  of  the  medium  such  as  variability  of  grain 
size  and  satmation  levels. 

Dick  presented  results  from  dynamic  loading  by  the  split-Hopkinson  pressme  bar  on  dry 
and  partially  saturated  Big  Black  and  Eglin  sands.  Results  were  that  dry  sand  crushes 
readily  and  void  collapse  was  incomplete  for  the  stresses  achieved  in  these  tests.  For  wet 
sands,  void  collapse  required  at  least  100  MPa  and  20  msec  before  significant  stresses 
could  be  transmitted.  Additional  measurements  of  pressure  in  wet  sand  as  a  function  of 
range  from  the  explosive  charge  indicated  that  the  stress  at  50  mm  range  (240  MPa)  was 
reduced  to  100  MPa  at  140  mm  range.  Chitty  and  Walker  presented  field,  laboratory, 
and  computational  investigations  of  wave  propagation  in  saturated  sands.  The  laboratory 
testing  focused  on  hydrostatic  compression,  uniaxial  strain,  and  triaxial  compression  tests 
to  measure  the  mechanical  properties  of  sands  to  support  their  effective  stress  modeling. 
Measmements  were  made  of  the  pore  fluid  flow  rates  through  various  sands  and  various 
means  to  measure  saturation  levels  in  situ  were  discussed.  They  concluded  that  the  mul¬ 
tiphase  effective  stress  model  reproduced  the  basic  physics  of  the  groundshock,  detailed 
skeleton  models  were  needed  for  accurate  late  time  simulations,  and  relative  flow  models 
were  needed  to  calculate  late  time  fluid  redistribution.  England  compared  measmred  water 
pressure  and  sand  stress  histories  for  a  scaled  1000  lb  burst  buried  25  ft  in  sand  under  50 
ft  of  water  with  results  calculated  with  an  effective  stress  hydrodynamic  model.  Predicted 
and  measmed  peak  water  pressures,  as  a  fimction  of  range,  agreed  to  within  30%  for  both 
fully  satmated  and  partially  satmated  sand  bottoms.  He  concluded  that  gas  filled  porosity 
is  the  dominant  material  parameter  controlling  shock  attenuation. 

The  micromechanical  models  were  presented  by  Pandy,  Shukla,  Barsoum,  and  Peric’. 
Pandy  presented  results  firom  a  percolation-based  computer  model  to  study  the  perme¬ 
ability  of  fliud  flow  through  porous  media  and  results  from  a  direct  simulation  Monte 
Carlo  model  to  study  shock  propagation  through  porous  media.  The  permeability  showed 
a  nonmonotonic  dependence  on  the  fluid  concentration  and  on  the  driving  bias.  The 
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shock  propagates  in  a  drift  like  fashion  under  high  porosity;  the  shock  front  damps  out 
very  quickly  with  low  porosity.  Shukla  discussed  experimental  and  numerical  studies  on 
stress  wave  propagation  through  small  polycarbanate  disks.  Results  showed  no  transverse 
propagation  of  the  stress  waves  through  linear  arrays  of  disks  under  load.  Propagation 
through  arrays  of  randomly  oriented  disks  of  random  size  was  dependent  on  the  structure 
of  the  material  fabric.  When  water  was  placed  between  the  disks,  the  magnitude  of  the 
contact  load  was  severely  attenuated  as  a  function  of  propagation  distance.  The  results 
were  compared  with  a  discrete  element  model.  Barsoum  discussed  a  microplane  approach 
that  included  strain  rate  effects  based  on  three  phase  micromechanical  behavior.  Peric’ 
discussed  the  importance  of  the  shear  band  instability  in  soil  mechanics  and  a  bifurcation 
in  the  pore  pressure  which  sets  limits  in  the  range  of  validity  for  the  expressions  for  the 
rate  of  displacement. 

The  synthesis  of  macroscale  constitutive  relations  from  microscale  constitutive  relations 
was  discussed  by  Flippen  and  Zhang.  Flippen  discussed  the  Condensation  Model  Reduc¬ 
tion  (CMR)  approach  which  can  be  directly  coupled  to  current  finite  element  models.  This 
model  has  achieved  some  success  in  modehng  laminates  and  ceramics.  He  discussed  imple¬ 
menting  a  macroscale  semi-discrete  constitutive  operator  synthesis  for  saturated  sediments. 
Zhang  employed  an  ensemble  phase  averaged  method  to  develop  constitutive  laws  for  FEM 
and  CFD  codes.  He  argued  that  this  method  provides  a  link  between  micromechanics  and 
the  macroscopic  model  while  retaining  the  physics  of  the  phase  interaction. 

V.  COMMENTS  ON  PRESENTATIONS 

It  is  clear  that,  to  lowest  order,  the  gas  content  is  the  dominant  parameter  controlling 
shock  attenuation  in  the  surf  zone  environment.  Even  a  minute  amount  of  air  drastically 
alters  the  pore  pressure  response  in  multiphase  porous  materials.  Air  content  as  small  as 
1%  can  reduce  the  shock  speed  by  factors  of  3  -  4,  and  reduce  the  shock  amplitude  by 
factors  of  20  or  more.  This  is  borne  out  by  both  experiment  and  numerical  simulation. 
It  is  also  clear  that  there  is  no  accepted  method  to  accurately,  and  rapidly,  measure  the 
gas  content  in  partially  saturated  sands/soils  either  in  situ  or  under  laboratory  controlled 
conditions. 

Several  of  the  computational  hydrodynamic  models  incorporated  continuum  models  for  the 
partially  satmated  sand  environment.  These  included  the  effective  stress  model,  the 
model,  and  the  multiphase  effective  stress  model.  The  predictions  from  the  effective  stress 
model  agreed  reasonably  well  with  experimental  results  (to  within  30%)  under  conditions 
of  relatively  large  loads  stemming  from  deeply  buried  mines.  This  is  to  be  expected  as 
the  controlling  parameter  with  large  amplitude  shock  waves  would  be  the  air  content,  and 
the  details  of  the  skeleton  structmre  of  the  sand  would  play  a  small  role  as  the  stresses 
are  well  beyond  the  yield  point  of  the  sand  skeleton.  The  Pq  model  reproduced  the 
experimentally  measmred  impulses  under  conditions  of  moderate  loads,  but  had  difficulty 
reproducing  the  peak  pressures.  The  details  of  the  shock  propagation  are  affected  by 
the  characteristics  of  the  sand  skeleton  which  are  not  modeled  with  the  Pq.  model  and 
the  Pa  model  does  not  include  strain  rate  effects.  The  multiphase  effective  stress  model 
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incorporated  compressibility  parameters  and  mechanical  properties  of  dry  sand  determined 
from  hydrostatic  compression,  uniaxial  strain,  and  triaxial  compression  tests,  which  lent 
some  validity  to  the  results.  The  numerical  results  for  partially  saturated  sand  looked 
physically  reasonable  but  there  was  no  comparison  with  experiment. 

These  continuum  models  appear  to  have  modeled  the  qualitative  features  correctly,  and 
the  quantitative  features  reasonably  well,  under  conditions  of  high  stress  and/or  high  air 
content.  It  is  unclear  whether  these  models  have  the  same  predictive  capability  under 
low/moderate  stress  and  saturated,  or  nearly  satmated,  conditions.  Under  these  circum¬ 
stances,  the  microscopic  structure  of  the  sand  skeleton  may  be  important.  The  microme¬ 
chanical  results,  both  experimental  and  numerical,  presented  at  the  workshop  indicate  that 
the  shock  propagation  characteristics  depend  strongly  on  the  microstructure.  These  results 
were  based  on  low/moderate  stress  levels  -  primarily  elastic  waves.  These  results  may  play 
an  important  role  in  sophisticated  explosive  net  designs  (multiple  shocks,  shaped  charges) 
and  could  play  a  role  in  the  long  time  evolution  of  a  single  shock  where  the  precursor  wave, 
caused  by  the  rigidity  of  the  sand  skeleton,  may  be  important. 

A  recmring  problem  facing  the  experimental  commtmity  is  the  accurate  measurement  of 
the  parameters,  and  the  controllable  repeatability  of  the  parameters  in  a  laboratory  setting, 
which  may  impact  the  evolution  of  the  shock  wave  :  pore  size,  grain  size,  density,  water 
content,  air  content,  packing,  compression  strength,  shear  strength,  anisotropy,  and  min¬ 
eralogy.  It  is  especially  difficult  to  control  the  degree  of  saturation;  although  the  limiting 
cases  of  dry  sand  and  fully  saturated  sand  (created  under  conditions  of  reduced  pressure) 
are  attainable.  It  is  nearly  impossible  to  produce  a  particular  saturation  level  a  priori  - 
again  pointing  out  the  need  for  a  rapid  and  accmate  means  of  measuring  air  content  in  wet 
sands.  Even  if  the  above  set  of  parameters  were  measurable,  to  a  reasonable  degree  of  accu¬ 
racy,  it  is  still  not  clear  what  role  the  above  parameters,  aside  from  air  content,  play  in  the 
evolution  of  the  shock  wave;  yet,  the  macroscopic  equivalents  of  these  parameters  are  key 
elements  to  most  computational  models  and  they  control  the  evolution  of  the  stress  waves, 
through  the  constitutive  relations  and  conservation  equations,  in  a  numerical  simulation. 
The  macroscopic  parameters  in  a  numerical  model  include  yield  strength,  shear  modu¬ 
lus,  plastic  strain,  work-hardening  coefficients,  longitudinal  and  transverse  sound  speeds, 
bulk  modulus,  Griineisen  constant,  and  Bauschinger  parameters  (controls  the  unloading 
characteristics).  These  parameters  are,  in  turn,  related  to  the  physical  variables  discussed 
above. 

The  split-Hopkinson  bar,  hydrostatic  compression,  uniaxial  strain,  and  triaxial  compression 
tests  are  good  first  steps  in  this  endeavor  to  measure  the  physical  characteristics  of  dry  and 
wet  sand.  Sample  preparation  for  these  tests  is  difficult  and  repeatability  is  questionable. 
The  Hopkinson  bar  test  is  based  on  one  dimensional  analysis  but  sand  grains  are  large 
relative  to  the  sample  size.  Early  evolution  of  the  stress-strain  relation  is  not  reliable,  yet 
this  region  is  important  because  the  air  voids  are  being  compressed  during  this  period. 
The  unloading  period  is  questionable  in  that  it  is  not  clear  whether  the  observed  hysteresis 
phenomena  is  due  to  plasticity  or  failure  (loss  of  water)  of  the  sample.  Crushing  of  the 
sand  sample  was  observed  for  the  dry  tests  but  not  for  the  wet  sand  tests. 
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The  integration  of  microscale  phenomena  into  macroscale  phenomena  amenable  to  a  con¬ 
tinuum  model  is  still  in  its  infancy,  although  homogenization  procedures  have  worked  quite 
well  for  ceramics  and  laminates.  The  reliability  of  these  procedures  rests  more  on  the  mi¬ 
croscale  model  chosen  than  the  procedure  itself  and  at  this  stage  it  is  not  apparent  which 
microscale  model  to  employ. 

VI.  RECOMMENDATIONS 

The  focus  of  this  workshop  was  on  the  constitutive  modeling  of  air  and  water  saturated 
sandy  media  as  it  pertains  to  the  development  of  accurate,  efficient,  physics  based  consti¬ 
tutive  models  for  gassy  sediments  in  and  near  surf  zone  environments.  The  constitutive 
models  control  the  evolution  of  the  shock  wave  and  improved  constitutive  models  are 
needed  if  progress  is  to  be  made  in  performance  prediction  of  mine  countermeasme  efforts. 
The  present  numerical  models  appear  to  model  the  qualitative  features  reasonably  well 
under  conditions  of  high  stress  and/or  high  air  content.  Detailed  comparisons  with  experi¬ 
ments  show  some  discrepencies  with  respect  to  peak  stresses,  shock  rise  times,  propagation 
speeds,  and  mine  damage.  Errors  and/or  inadequacies  in  the  constitutive  models  would 
manifest  themselves  in  this  manner,  under  the  assumption  that  the  conservation  equations 
are  being  solved  correctly. 

As  noted  above,  the  gas  content  is  the  dominant  parameter  controlling  shock  attenuation 
yet  there  is  no  accepted  method  to  accurately,  and  rapidly,  measure  the  gas  content  in 
partially  saturated  sands/soils  either  in  situ  or  under  laboratory  controlled  conditions.  It 
is  our  recommendation  that  ONR  take  the  appropriate  steps  to  address  this  issue  with 
the  following  two  caveats.  (1)  Numerical  simulations  indicate  that  little  additional  shock 
attenuation  occurs  when  the  gas  content  exceeds  5%.  If  the  distributed  explosive  net, 
presently  being  considered  for  mine  clearance  in  the  smf  zone,  can  be  designed  to  provide 
a  large  enough  pressme  load  to  trigger  or  disable  a  buried  mine  in  a  partially  saturated 
environment  with  5%  air,  then  there  is  probably  no  need  to  know  the  degree  of  satmation 
to  high  accuracy.  (2)  It  would  seem  possible  to  design  a  double  pulse  explosive  device 
with  the  intent  that  the  first  pulse  would  be  just  large  enough  (perhaps  in  the  1-2  kilobar 
range)  to  collapse  the  air  bubbles  in  the  medium  and  the  second  pulse  would  then  propagate 
nearly  rmattenuated  through  a  fully  saturated  sand- water  environment.  The  strength  of 
the  first  pulse  may  be  a  strong  function  of  the  gas  content  as  the  volume  strain  required 
for  pressure  saturation  of  the  pore  air  is  directly  related  to  the  initial  air  content  which 
varies  as  a  function  of  depth.  In  this  case,  knowledge  of  the  degree  of  saturation  would 
be  vital.  In  either  scenario,  the  key  engineering  goal  is  to  design  the  shape  of  the  shock 
wave  through  the  design  of  the  explosive.  The  sand  constituents  themselves  will  dictate 
the  shock  shape  required.  The  weak  link  here  is  knowledge  of  the  sand  constituents  and 
their  control  over  the  evolution  of  the  shock  wave. 

A  significant  effort,  entailing  significant  computer  resources,  should  proceed  to  define  the 
physical  variables,  and  the  parameter  range  for  those  variables,  that  influence  the  evolution 
of  stress  waves  in  a  wet  sandy  environment  under  low,  moderate,  and  high  loads.  This 
nmnerical  effort  should  be  correlated  with  an  experimental  program  based  on  stress-strain 
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response  measurements  and  shock  propagation  measurements  in  a  controlled  laboratory 
setting.  A  convenient  way  to  evaluate  the  constitutive  relations  is  to  do  experiments  and 
simulations  on  dry  and  fully  saturated  specimens.  This  data  could  be  used  in  estimating 
parameters  for  in  situ  materials  associated  with  a  given  event.  It  would  appear  to  be 
critical  to  accmately  model  the  stress-strain  response  of  the  ’’dry”  sand  skeleton.  The 
rigidity  of  the  skeleton  determines  the  velocity  of  the  precmsor  wave  and  the  elastic  limit. 
For  partially  saturated  sands,  the  material  strength  is  strongly  influenced  by  the  pore 
pressure  response,  which  is  controlled  by  the  volumetric  deformation  characteristics  of  the 
porous  skeleton. 

We  feel  there  is  a  need  for  a  more  focused  effort  in  micromechanical  modeling  coupled  with 
well  controlled  laboratory  experiments.  For  example,  propagation  of  stress  waves  through 
individual  grains  of  sand  is  probably  not  important;  quantiflcation  of  the  contact  forces 
between  grains  of  sand  is  probably  important.  If  it  appears  that  the  microstructure  plays  a 
strong  role  in  the  shock  evolution,  then  further  work  on  the  detailed  characterization  of  the 
surfzone  environment  is  required  in  addition  to  developing  continuum,  or  continuum-like, 
constitutive  relations  based  on  the  micromechanics  of  the  surf  zone. 

Since  the  ultimate  goal  is  a  viable  mine  countermeasure  program,  it  is  important  to  con¬ 
sider  extensions  of  the  research  effort  suggested  above.  It  is  our  recommendation  that  a 
focused  experimental,  both  in  situ  and  laboratory,  and  numerical  simulation  program  be 
developed  to  investigate  both  the  single  pulse  and  multiple  pulse  explosive  net  designs  in 
more  detail  with  the  goals  being  (1),  a  quantitative  description  of  the  shock  propagation 
characteristics  as  a  function  of  gas  content  and  explosive  load  for  both  explosive  net  de¬ 
signs  and  (2),  a  prediction  of  performance  capability  for  both  distributed  explosive  net 
designs.  The  following  issues  need  to  be  addressed  :  experimental  measmement  of  shock 
attenuation  for  air  content  exceeding  5%  and  comparison  with  computational  results;  ex¬ 
perimental  data  for  shock  propagation  characteristics  in  a  fuUy  saturated  sand  environment 
as  a  function  of  range  and  comparison  with  computational  results;  numerical  simulations 
detailing  the  strength  and  dmation  of  the  loads  required  for  pressure  satmration  of  the  pore 
air  as  a  function  of  air  content  and  range  with  experimental  veriflcation;  and,  numerical 
simulations  of  the  loads  and  timings  required  for  the  double  pulse  explosive  net  design 
with  experimental  veriflcation. 

In  summary,  our  recommendations  to  NRL  and  ONR  are  : 

•  make  a  concerted  effort  to  develop  techniques  (both  laboratory  and  in  situ) 

to  accurately  and  rapidly  measure  air  content  in  partially  saturated  sandy  media, 

•  develop  a  detailed  experimental  and  computational  program  to  deflne,  measure, 
and  calculate  the  physical  variables  controlling  shock  evolution  in  a  wet, 
sandy  environment, 

•  determine  which  micromechanical  properties  are  important  to  a  physics  based 
constitutive  model  and  develop  the  means  to  incorporate  those  properties  into 
continuum  numerical  (hydrodynamic)  models. 
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•  enlarge  the  scope  of  distributed  explosive  net  designs  to  include  multiple  shocks 
and  shaped  charges  with  computational  and  experimental  support. 

Several  of  these  efforts  may  already  have  taken  place  or  are  now  taking  place  within 
the  ONR  community,  but  this  was  not  apparent  at  the  workshop. 
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